Temperature shifts in the Sinai model: static and dynamical effects 
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We study analytically and numerically the role of temperature shifts in the simplest model where 
the energy landscape is explicitely hierarchical, namely the Sinai model. This model has both 
attractive features (there are valleys within valleys in a strict self similar sense), but also one 
important drawback: there is no phase transition so that the model is, in the large size limit, 
effectively at zero temperature. We compute various static chaos indicators, that are found to be 
trivial in the large size limit, but exhibit interesting features for finite sizes. Correspondingly, for 
finite times, some interesting rejuvenation effects, related to the self similar nature of the potential, 
are observed. Still, the separation of time scales/length scales with temperatures in this model is 
much weaker that in experimental spin-glasses. 



I. INTRODUCTION 

The phenomenology of very different glassy systems (spin glasses, structural glasses, pinned 
defects) exhibits striking similarities. This is often discussed in terms of complex energy 
landscape, with some universal features in the organisation of valleys, saddles and barriers 
[0, H, H, § . Mean-field models suggest that the energy landscape of some systems are hier- 
archical, with a self-similar organisation of valleys within valleys Q . This picture has been 
advocated to interpret 0, ^ the striking rejuvenation and memory effects during temper- 
ature cycling experiments, observed first in spin-glasses, but also in disordered ferromagnets, 
dipolar glasses, polymer glasses (PMMA) or gelatin ||, [l^, 11, 



Another general framework to understand disordered systems is the droplet picture, de- 
veloped in the context of spin-glasses and of pinned manifolds (domain walls, vortices, etc.) 
[ 1^, . The approach aims at describing the physics in terms of localized excitations of 
different size (droplets), with an energy scale that grows with their size. As emphasized in 
[ 15| , [l^ , these two approaches are in fact complementary, and the dynamics in the droplet 
model is naturally hierarchical, due to the strong separation of time scales with length scales 
and with temperature. Temperature plays the role of a microscope since the active length 
scales that contribute to the dynamics are very different at different temperatures. 

In both pictures, memory conservation is due to this separation of time scales with tem- 
perature 0, |l| (see also |l|]). The origin of rejuvenation, however, is different. In the 
hierarchical landscape framework, a small temperature drop (from Ti to T2 = Ti — AT) 
reveals finer details of the random energy landscape, among which the system must reequi- 
librate. In other words, small length scales, that are 'unpinned' at temperature Ti, freeze at 
T2, thereby producing a strong out of equilibrium signal (rejuvenation). This scenario can 
be given some precise meaning in the context of the Random Energy Model (rem) close to 
its critical temperature pC[ |, or in the generalized (multilevel) rem where rejuvenation and 
memory effects very similar to experiments can indeed be observed pl| , p2| . On the other 
hand, rejuvenation effects in the droplet model have been attributed to 'temperature chaos': 
for any temperature difference AT, the equilibrium configurations at the two temperatures 
are completely different beyond a certain length scale ^at, called the overlap length, that 
diverges when AT ^ as a power-law. Correspondingly, the correlation of the free-energies 
at the two temperatures are thought to decay with the size of the system L as exp(— L/^at)- 



This scenario, postulated for spin-glasses in [|13|, 14 , has been given credit recently in the 
context of directed polymers psf |2J, In this case, both analytical arguments and nu- 
merical results point towards the existence of an overlap length, although the free-energy 
decorrelation appears to be much slower than exponential in the size of the system. For 
spin-glasses, the status of temperature chaos is still controversial and unclear as compared 
to the influence of other, stronger perturbations such as magnetic field changes or coupling 



2 



strength changes. Results in mean-field spin glasses ^ and numerical calculations for 
small systems hint at extremely weak chaotic effects in temperature. This conclusion 
extends also to short-range systems where chaotic effects in the equilibrium properties seem 
to be extremely weak [ pOj |30f| . Chaotic effects, if any, seem to appear only for quite large 
systems even when AT/T is of order one Nevertheless, some recent spin-glass exper- 

iments were interpreted in terms of an overlap length [^3| , whereas numerical simulations of 
the 4D Edwards- Anderson model has revealed rejuvenation due to small scale freezing with- 
out chaos . (For a detailed study in the 3D Edwards- Anderson model and a discussion 
about why rejuvenation effects are not observed in simulations see Q). 

In view of the controversy, we feel that more work on the subject is needed. The aim of the 
present paper is to study in details the question of 'temperature chaos' on the Sinai model, 
which is the simplest model where the energy landscape is explicitely hierarchical, and for 
which many exact results are available H, |36[ |3^, ||, ^ . The Sinai model was discussed in 
this context in |l^; preliminary results were obtained by H. Yoshino, but no systematic 
study had been performed. 

The Sinai model is an example of a one dimensional self-similar potential with long-range 
correlations. In high dimensions, this problem is equivalent to a mean-field spin-glass with a 
continuous Replica Symmetry Breaking solution In one dimension, however, there 

is no phase transition: the long time, large scale behaviour of the system is ruled by the 
zero-temperature fixed point, where the deepest minimum dominates. This is because the 
entropy in this model is of order one, whereas the energy scale grows as ^/L. Since the 
location and energy of the minimum is temperature independent, there cannot be any true 
temperature chaos in this problem. However, there are interesting static and dynamical 
crossover effects, as a function of the size of the system or of the time, that we study in 
details in this paper. We find that rejuvenation and memory effects are present, in embryo, 
in the Sinai model. 

The paper is organized as follows. In the next section we present the Sinai model. In 
section III we analyze the statics of the model. In section IV we discuss the dynamics, first 
by studying aging in the correlation length and then by studying the a.c. susceptibility 
during temperature cycles. In section V we give the conclusion derived from our analysis. 



II. THE MODEL: SINAI POTENTIAL 



The Sinai model belongs to a wider class of random potential models. It describes the 
dynamics of a point particle under the action of a random, uncorrelated force field F{x) 
which models several physical situations ]3^ , such as the motion of domain wall in the 
Random Field Ising Model or the motion of a dislocation kink. More recently, this 
model was argued to be relevant to describe some tapping experiments in sandpiles p3[ 
and the unzipping transition of DNA Q]. Interestingly, the effective potential acting on 
the end point of the disordered directed polymer in 1-1-1 dimensions is also of the Sinai 
, 0. However, in this case, the effective potential itself becomes temperature 
dependent, and the role of temperature changes in the directed polymer is much more subtle 

II, n, n. 

In the following, we will study the discrete Sinai model where each site corresponds to a 
different state or configuration with energy Vi. The system consists of a box of length L in 
which we generate a random potential. Each sample of the random potential is constructed 
as follows: at each site i e [1,^], we generate a random Gaussian force fi with zero mean 
and variance fifj = Sij. All along this paper, we have taken tr = 1 in the numerical 
simulations. The potential in each site is the sum of the forces in the previous sites Vi — 
— '^j^ii Fj, and is thus a random walk as a function of the position. Thus correlations 
(and equivalently barriers) grow like: 

W^W = ^'\^-j\- (1) 

This model has no thermodynamic transition so that in the limit L oo, the physics is 
dominated by the T = glassy fixed point |^, |4^. As is obvious from Eq.(|l]), a rescahng 
of the length by a factor b is equivalent to a change of the scale of the potential by a factor 



type m 46 



3 



or of the temperature by a factor 1 / \/h. As far as the statics are concerned, this means 
that being at low temperature in a small system is equivalent to having a larger system but 
a higher temperature. For thedynamics a change in temperature leads to both a change in 
length scales and timescales. 



III. STATIC CHAOS 



A. Observables 



Our aim is to investigate how the thermodynamic properties of the Sinai model change 
when comparing the same system at two different temperatures. In order to probe the change 
in the free-energy landscape we have studied two different quantities, (i) the correlation 
function of the free-energy fluctuations at different temperatures and (ii) the correlation of 
the particle position. More precisely, we have studied the following observables: 

• Free-energy correlations: It consists in measuring the free energy fluctuations of the 
system at different temperatures averaged over the disorder: 



where Ft is the free energy at temperature T and (..) stands for the average over 
different realizations of the potential. This quantity was originally proposed by Fisher 
and Huse |^ to study temperature chaos in the directed polymer problem. The fact 
that this correlation tends to zero with the size of the system shows that different energy 
valleys contribute to the total free-energy at different temperatures. More recently this 
quantity has been used to study the chaotic properties in the REM [Q, and in the 
directed polymer problem when different type of perturbations are introduced Q , or 
for directed polymers on a hierarchical lattice |^ . 

Fluctuations of position: A more geometrical way to visualize 'temperature chaos' 
which could be of direct interest in some cases, for example the zipping and unzipping 



problem of DNA |44 , is to consider the following quantity. To each site i, we associate 



a position Xi = i/L and compute: 



di,2 ^ ((x)t, - (x)tJ2, (3) 

where (..)t is the thermal average at temperature T, and (...) is the disorder average. 

The study of the distance between average positions corresponding to different tem- 
peratures gives us an indication of the distance between states contributing to the 
partition function Z at different temperatures. If typical states contributing to 2, at 
Ti and T2 are completely different, then (ii,2 will remain finite as L ^ 00. Note that 
di_2 has an upper bound d\;2. < 1/6, where 1/6 is the value reached if the occupied 
sites are completely uncorrelated. 

Both quantities Cp and d have been studied in the rem to show that even if the 
energy landscape is fixed (in the sense that there is no reshuffling associated with the different 
valleys as temperature is changed) , the model exhibits temperature chaos around the critical 
temperature, where the Boltzmann weight 'condensates' into a finite number of sites. Because 
there is no finite temperature transition in the Sinai model, one only observes mild effects 
under a temperature change, that are maximum around the crossover temperature T ^ L^/^ 
(see below). 

Numerically we have investigated how these quantities behaves as a function of L for 
different temperature changes. We have averaged over 2000 potential samples of sizes ranging 
from L = 1 to L = 2^^. 
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B. Thermodynamics of the model 



The thermodynamics of the Sinai model has been well studied |48, In the large 

L limit temperature is irrelevant. The system is frozen or localized in the minima of the 
potential, so that physical observables are governed by the ground state and its fluctuations. 
The free-energy grows as ~ — V^i independently of temperature and the entropy reaches 
an L-independent value which depends on temperature S = A + 2 In T. The free-energy 
fluctuations are also essentially dominated by the fluctuations of the deepest valley at 
any temperature. Therefore, it is reasonable to expect that in the thermodynamic limit a 
change in temperature has no significant effect on physical observables since the ground 
state properties are temperature independent. 

The finite size corrections to the thermodynamic behaviour enter through the variable 
g — (7^/L/T. Corrections are important when g ~ 1, or L ~ L* = {T/a)'^, signaling the 
crossover between two different limits: 

• L ^ L*: energy differences are much smaller than the temperature, so that all the 
sites contribute to the partition function. 

• L ^ L*: energy differences and barriers are huge and the particle is localized; only a 
few sites within the deepest valley contribute significantly to the partition function. 

This crossover is clearly observed in the behavior of the fluctuations of the mean position 
(in units of L): 
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= (a;)2 - (x) . (4) 

In figure |^ we show the results for different temperatures versus the scaling variable g 
defined above. The small g behavior can be computed using a high-temperature expansion 
of the partition function: 

A' = ^+0{g^). (5) 

As we can see in the figure, the behavior for small g matches nicely this prediction. For 
large g the system is completely governed by the T — behavior, thus average position 
fluctuations should approach the fluctuations of the ground state position. In figure ^ we 
also plot the ground-state position fluctuations as a function of VX, that converges towards 
1/12. Numerically we can see that this crossover takes place at g ~ 3, corresponding to 
L* ~ 10(T/cr)2 

This static crossover can be mapped to the dynamical crossover from a Brownian diffusive 
regim e to an activated regime as we will see in the analysis of the correlation length in section 



IVA. 



1. Free-energy fluctuations 

As we can see from the inset in figure ^, the free-energy correlations are a non-monotonous 
function of L . For small L the free-energy at different temperatures start decorrelating 
but at a certain length, Lp(Ti,T2), Cp reaches a minimum. Then, it starts increasing back 
and reaches Cp = 1 when L oo regardless of the temperature difference (provided it is 
finite) as expected. One can distinguish between two different regimes: 

For small system sizes/large enough temperatures the energy landscape is essentially 
flat, so that all the sites contribute to Z. Therefore, in this regime a high temperature 
expansion is expected to yield the correct behavior of Cp . The result is the following: 

1-Cf cx {(3i- f32f<7^L (A/3crVI^0), (6) 



where P=l/T. 
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FIG. 1: Sinai Potential: plotted versus the scaling variable g = a\fLjT for T = 

1, 2, 10, 50 and 100. The dashed line corresponds to the position fluctuations at zero temper- 
ature. For small g we have plotted the high temperature prediction 3^/120. 



• L > LJ, 

In this limit the system is governed by the ground state and its fluctuations, thus we 
expect that for L ^ oo the correlation is perfect. However when gi^2 = o"\/i/Ti^2 is 
large but finite, the free energy decorrelates slightly. The behavior in the large L limit 
can be understood by very simple arguments. Suppose that T = (Ti + T2)/2 is small, 
and that the relative temperature change is also small: e = (Ti — T2)/T ^ 1. Using 
dF/dT = ~S{T), where S{T) is the entropy at temperature T, one finds: 

AF = F{T2) - F(Ti) = AT S{T) + 0{e^), AT = Ti - Tj. (7) 

Note that this relation is also true for the fluctuating part $ of F. Substituting 
into the expression for Cf and expanding for small e, we finally get: 

Cf^1-AT^^^^.. (8) 

We show in Appendix A that the entropy fluctuations tend to a numerical constant K 
that we compute in the large L in limit. Thus recalling that free-energy fluctuations 
scale as $2 = K'a^L, where K' can be computed from [Q, we finally find the following 
scaling behavior for 1 — Cf'- 

KAT^ 

• Crossover 

The crossover between both regimes will take place at a certain length Lp such that 
APay^ - AT/o-yZ*^, which yields L*p - TiTa/cr^. The maximum of 1 - Ci^ 
actually occurs for Lp « 10 TiT2/a'^. This suggests the following scaling form for 
I-Cf- 
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with f{x) ~ X for a; ^ and f{x) 1/x for a; — > oo. The resulting scahng plot is 
shown in figure ^. This scaling is acceptable only when temperature differences are not 
too large, as one would expect. 



(1-C,) l,/at' 




FIG. 2: Scaling function for the free-energy correlations. (1 — Cf) x \j,2 as a function of -pr for 
different pairs of temperatures {Ti,T2). 



2. Position fluctuations 

In figure ^ we show the curves corresponding to the average distance shift di 2 defined by 
Eq.(||) above, for several temperature differences. Again we find that this quantity is a non 
monotonic function of L. As a matter of fact, the behavior of both di,2 and Ci? is similar, 
and the same type of argument can be used for both quantities. 

• L<L2 

In figure ^ we can see that for large temperatures, di,2 increases with system size. This 
behavior can be explained by looking at a high temperature expansion. As shown in 
Appendix B, we find in this regime di^2 ~ ■j^i^f^^'^^L. 

In the large L region the behavior of di,2 should be governed by the ground state 
fluctuations. One starts from the relation: 




where (...)t is the thermal average at temperature T. This relation is derived from 
the fluctuation-dissipation relation {d{x)/d(3 — ~{xVx)c) and integrating it assuming 
e <C 1. The main contribution to c?i_2 for small temperatures will come from two nearly 
degenerate valleys that are a distance ~ L apart. Restricting the thermal averages to 
these two valleys, one finds: 

where D is the distance between the two valleys, and (jrjy/L is their energy difference, 
such that ry is a positive random variable of order unity. In order to obtain di^2 we have 
to compute {{X2) — (xi))^. This average depends on the joint probability distribution 
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of the excitations P{ri, D). It is reasonable to assume that the probabiHty distribution 
factorizes as P{r], D) = h{r]) p{D). Therefore in the Hmit /3a\/L ^ oo we have: 



/•= 

Jq 



drj h{ri)e 



-2f3<T' 



n 



(13) 



RecaUing that x is the rescaled distance i/L, one expects D"^ ^ 1. Provided h{Q) ^ 0, 
one finally finds: 



di_2 ~ ^ {(3aVL oo). 



(14) 



• Crossover 

As for the free-energy correlation, we can extract a crossover length scale which sepa- 
rates both regimes. The crossover length that is obtained is ^ (Ti T2)^T~3 ^ that 
coincides with Lp in the limit Ti — T2. Again, we can try a scaling formula: 



(AT)^ 
(TiT2T)i 



where 



7(x)^ 



■' X 

1 



X 

X - 



^ 

00 



(15) 



Note that the functions / and g are quite different. This reflects the fact that that the 
two observables probe different mechanisms. In figure || we show the scaling plot for 
different pairs (ri,T2). Notice that all the curves display the maximum at the same 
value of the scaling variable -pr — 10. 



C^: T,=10 V50 
□ — □ T,=1 .5 V5 
T,=0.5 V2 
T,=1 V2 
T,=1 V 5 
T,=1.5 V10 

□ — qTj.s V 10 




FIG. 3: Rescaled average distance shift di,2 x ''^(^^rp ^ plotted versus L/L*^ for different T\ and T2. 



C. Discussion 



The outcome of the numerical analysis of the Sinai potential is clear: in the thermodynamic 
limit statistical properties are governed by the T = fixed point. This means that regardless 
of temperature, for large enough system sizes, the system only sees the free-energy valley 
associated to the global minimum. Effectively this situation is equivalent to say that in 
this limit there is no chaos in temperature since the statistical properties are those of the 
minimum of the potential. 
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This situation is very different from what happens in the directed polymer problem. This 
model, as well as the Sinai model, has no thermodynamic transition. It only displays a 
low-temperature glassy phase. In a recent paper this model has been shown to be 
extremely sensitive to any perturbation leadingto the vanishing of correlations between 
systems at different temperatures (see also [^s], |3^ for a related discussion.) In this model 
this is due to existence of anomalous large excitations which have a very low free-energy cost. 
These excitations cost a lot of energy AE ~ L^/^ but are very favored entropically, since 
T AS" ~ L^/^, so that these two contributions may cancel to yield a low cost in free-energy. 
In the Sinai model these anomalous excitations do not exist because the entropy is very small 
(not extensive) and never cancels the energy cost (~ L^/^) of such excitations. 

It is interesting to compare the above crossover lengths p with the crossover length at 
each temperature ^ (T/cr)^. The maximum decorrelation takes place for a system size 
such that L2 < L = L*^ p < L\. In this case the system at T2 is already localized whereas 
the system at Ti is still delocalized. The strong influence of temperature shift in this case is 
a smeared out version of the inflnite susceptibility found in the Random Energy Model |po| . 
In this case there is a true finite temperature phase transition, and not a mere crossover as 
in the Sinai case. For larger system sizes, when L\ < L2 < L, both systems are governed by 
the zero temperature fixed point and correlations increase. 



IV. DYNAMICS 

The dynamics of this model have been well studied, both analytically |^ , Q and numer- 
ically [ p7| . Single-time as well as two-time quantities have been analyzed. Here, we study 
one observable of both types; we define in particular an 'a.c. susceptibility' that should be 
closely related to the analogous observable studied in spin-glasses. 

In our simulations we have used boxes of length L = 1024 with periodic boundary con- 
ditions. The dynamics has been simulated by the Monte Carlo method using Metropolis 
algorithm. For each realisation of the random potential, in order to sample adequately the 
energy landscape we have considered all possible initial conditions (L = 1024) and we have 
averaged over n ~ 100 different histories for each starting point. The total number of sam- 
ples used in temperature cycling experiments is around 200 — 300. We have analyzed two 
different quantities: 



the correlation (or explored) length: 



(16) 



The brackets and overline mean that we average over both the L x n histories and 
samples respectively. The initial condition is a uniform distribution equivalent to 
a quench from infinite temperatures. This correlation length only gives information 
about the 'large scale' mechanisms and the temperature cycling effects on this quantity 
can be fully explained in terms of effective times. 

• In order to probe 'smaller' length scales that are more sensitive to temperature changes, 
we have defined the following 'a.c' susceptibility xi'-^itw)- 



(17) 



where the average is taken over the probability V{x,tw) that a particle is at position 
X at time t.^, with a uniform distribution of particles at time t — |59t] . In other 
words, we measure the typical extra distance travelled by particles during a time l/w, 
weighted by the dynamical distribution at time tw ■ Such a quantity was also considered 
in i^,|3^. 



Our main interest in this section is to study the temperature cycling experiments which 
are carried out in spin-glasses that have shown striking rejuvenation and memory effects 
Esl 16 . Our main goal here is see to what extent these effects are already present in the 
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Sinai model. We have performed numerically the standard temperature cycling experiment: 
quench from infinite temperature down to Ti and let the system relax during t^-^ ; then 
change the temperature to T2 =Ti + AT and let the system evolve during t^j^ ^.nd finally go 
back to Ti . We have studied cycles with positive and negative AT for several waiting times 
and frequencies. 



A. The correlation length 



The time evolution for the correlation length at different temperatures is shown in the left 
panel in figure IJ. The growth of the correlation length depends exclusively on a temperature 
dependent microscopic timescale to(TY This timescale is related to the crossover between 
two different dynamical regimes p5|, ^ : 



• t <ti tq: Short time dynamics where no barriers are present, so that in this regime we 
have usual Brownian diffusion £,'^{t); = D t. 

• t ^ tq: Long time, activated dynamics, with an activation time which follows an Ar- 



rhenius law with a typical barrier _B ^ (T\/L | 35 , leading to ^^(t) ~ (T/cr)^ ln^(— ) 



The crossover takes place when barriers become comparable to temperature so that activation 
between valleys dominates the dynamics. This crossover is directly related to the static 
crossover found for position or free energy correlations from a high temp erature regime (no 
barriers) to a (thermodynamic) low temperature regime (see section [II B). The microscopic 
timescale can be thus identified with the time that typically the system takes to explore 
this static crossover length scale L* ~ 100 (T/ct)^. Therefore, for cr = 1, to(T) = L*'^/D ~ 
200 T", (D can be evaluated to be £) « 0.5). 
One therefore expects that the correlation length at a given temperature ( p^ can be 





Fit with log correction 
Fit wittiDut log corrections 



FIG. 4: Left. Correlation length versus time for different temperatures. The solid line corresponds 
to the short time diffusive behavior ^(t) ~ \/?. In the inset we plot ^^/T* versus the scaling variable 
T = t/200 T*. Averages over 75 x 1024 histories. Right. \nt versus for T = 0.7 fitted with the 
effective function in (|l9|), with and without the lnln£ correction. 



expressed as T^/(t), where / is a function of the scaling variable r — t/ro, with: 



T « 1 /(r) 
T»ro(r) /(r)^ 



InV 



(18) 



This scaling behavior works very well, as shown in the inset of Figure ^ where we rescale 
together all temperatures. Such a crossover, between a short time growth law and a long time 
activated behaviour is also expected in the droplet description of spin-glasses ||l^, |l6, 18, [sH , 
or directed polymers, where barriers grow with the size i of the excitations as BJT) = TP , 
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where T is a function of temperature. This leads to a logarithmic growth of the size of the 
droplets, I ~ (In—)'?, where tq is a microscopic attempt time, possibly renormalized by 

critical fluctuations ||l6|, |5^, The Sinai case corresponds to -0 = 1/2 and tq oc As a 
more precise description of the crossover, we have shown in Fig. ^ the following fit: 



t{t) = Ae exp{B Vilnlni), (19) 

where the In In £ accounts for the famous Khinchin iterated logarithm law for the maximum 
of a random walk [Q. Note that the two limiting behaviors for in ( p!8| ) consistent with 
this last expression. In figure ^ the right plot shows \nt versus for T = 0.7 which are 
nicely fitted by the expression above. 

The correlation length is a monotonic function of time and temperature. The temperature 
only plays the role of slowing down dynamics, therefore a change in temperature only changes 
the growth law. Any time spent at a temperature Ti is equivalent to having spent an 
effective time tcff (^w, ^i, T2) at T2 such that ^^(t^j, Ji) = (icff , ^2 ) • Thus in temperature 
cycling protocols there is no trace of the chaotic effects observed experimentally on the 
correlation length itself, see Fig. 0. 
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T,=1^T2=5 
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FIG. 5: Correlation length ^(tm) versus time for the following experiments: ie/t) spend tw at 
Ti = 5 and then quench the system down to T2 = .5 for — 16, 1662, 8192, 40959 and 
106494 from bottom to top; Right) spend tw at Ti = 1 and then heat the system up to r2 = 5 for 
U = 79, 1662, 8192 and 106494. 



B. Susceptibility 

In the previous section we have studied the correlation length ^(^u,) that tells us how far 
the particle can go in a time tw This quantity is the analogue of the size of the domains in 
a droplet coarsening description. However, in a.c. experiments after a negative temperature 
shift one probably observes how the 'domain walls' rcconform on a scale which is small 
compared to Citw)- The susceptibility defined in ( p7[ ) probes these 'small' length scales 
and might show an interesting behaviour during temperature cycling, not revealed by ^{tw) 
(see Fig ^). The study of the 'response' function is useful because the results admit 
an intuitive interpretation in terms of the evolution of V{x,tw), which is the quantity that 
keeps track of the thermal history of the system. 

In order to compare with experiments, one should be in the following conditions: long 
waiting times and low frequencies u (as compared to microscopic timescales), but such 
that uj tw 3> 1. This last condition is imposed by the fact that a harmonic response can 
only be measured on a time larger than one oscillation period. This also ensures that one is 
in a regime where the violations of the Fluctuation Dissipation Theorem are weak and one 
can identify the fluctuation that we measure, (pT[), to a response [pT|. 
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FIG. 6: Measure of the effective time when coohng the system down to T2 — 0.5 after having spent 
twi ~ 512 at Ti = 5 at Lj = 1/128. Note that only the late part of the curves can be superimposed: 
there is a transient that cannot be accounted for using an effective time. 
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FIG. 7: Normalized effective times foffAmi: plot A) at fixed lu versus t^i, plot B) at fixed t^i 
versus l/uj. 



From the results of the simulations we observe that the effect of aging at temperature Ti 
on the relaxation at T2 depends strongly on the temperature difference and on the waiting 
time. This effect can be quantified by defining an effective time. For instance, when cooling 
the system from Ti to T2 one expects that if the system is completely rejuvenated, the 
relaxation curve XTi,T2 (i^j ^id2) should correspond to that obtained after quenching from high 
temperature Xoo,T2 (f^i ^«;2)- (Here t^2 is counted from the time at which the system reaches 
T2). However, if the relaxation at Ti affects aging at T2 then rejuvenation is only partial and 
the new relaxation corresponds to that of the system after aging during an effective time 
toff at T2, Xoo,T2 (1^: ^11)2 + tos). Thus if teff = 0, rejuvenation is complete. In Figure |^ we 
show how this effective time is measured. Note that only the late part of the curves can be 
superimposed: there is a transient that cannot be accounted for using an effective time. A 
similar effect can be observed in spin-glasses. The same effective time can also be defined 
when heating back the system, as a measure of memory recovery. 
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Experimentally, the observed facts in spin glasses are the following. 

• When cooling the sample from Ti < Tc to T2 < 7i, t^s is larger than t^i when AT is 
small enough, indicating that aging at Ti is more efficient because it is faster at higher 
temperature ||55|] . However, as AT increases, icff starts decreasing towards zero. In 
other words, for large enough AT, rejuvenation is complete. A related phenomenon is 
the absence of cooling rate effects. The relaxation at a low temperature does not depend 
on the cooling history. Only the very last steps at nearby temperatures determine the 
final relaxation [pl]| . 

• When heating back from T2 to Ti the system resumes its relaxation, for large enough 
AT, as if the stay at lower temperatures did not take place (perfect memory). As 
AT is reduced, an effective age < that accounts for the time spent at T2 must be 
added. Furthermore, a small transient appears at short times, that was called 'memory 
anomaly' in [55[| . Surprisingly, this anomaly is non monotonous with AT: for small AT, 
the reference curve corresponding to perfect memory is reached from below, whereas 
for larger AT, it is reached from above Finally, for still larger AT's, memory is 
perfect, as stated above. 

Note that if the time t^i spent at Ti in the first stage of the cycle is very small, one 
expects, from the previous discussion, to see 'rejuvenation' after a positive temperature 
shift, since the system has kept the memory of the age it had on the way down, i.e. 
twi <C tw2- This rejuvenation after a positive AT has been observed many times 
experimentally. 

Let us now turn to the results of the simulations, following the above presentation of 
experimental data. 

• The data of Fig ^ corresponds to temperature cycles with Ti = 5, T2 = 0.5 or Ti = 2, 
T2 = 0.5. The main observation here is that clear rejuvenation is observed, with an 
effective shift time that decreases as AT increases, as in the experiments. Intuitively, 
this corresponds to the fact that since the potential is self-similar, the local dynamics 
probed by x(w; ^u;) is not sensitive to the depth of the potential valley that is currently 
occupied. Therefore, aging at Ti has already selected some low-lying valleys, but the 
intra-valley dynamics is insensitive to this. This argument for rejuvenation based on a 
hierarchical energy landscape, and in the absence of temperature chaos has been put 
forward i n 0, |l^ |6j, and has been confirmed numerically in the multi-level trap 
model in | |21| . Lower frequencies, that correspond to larger length scales, are less easily 
rejuvenated, as expected, since the separation of time scales is not as sharp: see Figure 
0, where we show the normalized effective waiting times as a function of both t^i and 

• When heating back the system to the initial temperature Ti, some memory is observed. 
However (i) when the temperature difference is not very large, some effective time, 
accounting for the period spent at T2, must be included, as in the experiments, and 
(ii) a strong transient 'memory anomaly' is observed, even for quite large AT's (see 
plots A and B in figure Note that we have always observed this memory anomaly 
to be negative, i.e; the reference curve is reached from below. This memory anomaly 
is defined as: 

/^X^x{twi+tu,2 + l/^)-x{tZ,i) . (20) 

where xi^w = t~i) is the susceptibility just before the quench and xitwi + tw2 + 
corresponds to the first possible measurement at frequency lo after heating back to Ti . 
In Figs. H and ^ the dependence of the memory anomaly on different parameters, 
including the frequency, is shown. 

In fig. 1^ plot we show how the memory anomaly varies with the time tuj2 spent at 
T2, for a fixed uj = 1/128 and twi = 1024, for the largest AT. From the inset, we 
see that the larger AT, the smaller |Ax|, which agrees with the interpretation that 
the strong rejuvenation effect found for large AT arises from the separation of length 
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FIG. 8; Negative and positive temperature cycles. Plots A and B cycles with AT < 0, data for 
uj = 1/128, = 128, = 1024 and for cycle A: Ti = 5, T2 = 0.5 and B: Ti ^ 2, T2 ^ 0.5. Plots 
C and D cycles with Ar'> 0, for cycle C: Ti = 0.5, T2 ^ uj ^ 1/1024 and = t„2 = 8192 and 
D: Ti = 0.5, T2 = 5,Lu = 1/1024 and U,, = t„2 = 8192 



scales. On the plot in fig. we show how the memory anomaly depends on the 
frequency uj. Since the susceptibility at t^i itself depends on lu we plot the relative 
variation of the susceptibility with respect to xit^i) fixed t^i — = 1024. In 
the inset we show Ax/x(i^i) versus lj. Note that Ax/xi^wi) is always negative and 
decreases with increasing 1/uj. In p5| , Sasaki et al. find that in the anomaly can be 
both positive and negative when one works in the vicinity of a transition temperature. 
In the Sinai model we have not been able to observe such a positive anomaly. For 
smaller AT's these effects are blurred because length scale separation becomes weak. 
Aging at different temperatures is cumulative and Ax/x is also larger (see fig ^ for the 
Ti — 2 T2 — 0.5 cycle). For these smaller AT's, we have found that the memory 
anomaly becomes non monotonous with frequency. 

• In the positive temperature cycle experiments 

We have also observed that the third stage is completely independent of i^ji- Heating 
back the system to T2 erases the initial aging accumulated at Ti . This is expected, since 
the dynamics at T2 allows the system to leave the traps that it had slowly explored at 
Ti. The new relaxation at Ti is aged, but the age is only due to the effect of aging at 
T2. This is similar to the effect observed experimentally. 



V. DISCUSSION 



In this paper we studied in detail the role of temperature shifts in the simplest model 
where the energy landscape is explicitely hierarchical, namely the Sinai model. This model 
has both attractive features (there are valleys within valleys in a strict self-similar sense) , but 
also one important drawback: there is no phase transition so that the model, in the large- 
size limit, is effectively at zero temperature. Therefore, in this limit, temperature shifts do 
not lead to interesting phenomena: entropy in this model does not play any role, so that 
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FIG. 9: Memory anomaly in a negative temperature cycle at fixed ui. We plot the x — x(^™i) 
fixed Lu = 1/128 and t^l = 1024 for the temperature changes Ti = 5 = 0.1,0.2 for different 
tw2 ~ 512, 1024, 4096. In the inset we plot Ax as defined in (EOl) for the different tw2- 
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FIG. 10: Memory anomaly in a negative temperature cycle at fixed AT. Ax/xi'^wi) cycle 
Ti = 5 r2 = 0.2 with U„i = = 1024 we plot Ax/xi^Zi) for different frequencies, lo = 1/4.. ..1/256. 
In the inset we plot Ax/x versus l/co for this case and for the cycle Ti = 5 r2 = 0.1 with the same 
parameters. 



excitations have an enormous free-energy cost and cannot be favored by the temperature 
perturbation, contrarily to what happens in a closely related model, the Directed Polymer 
H (see also g). 

Nevertheless, for finite sizes/finite times, some interesting crossover phenomena qualita- 
tively reproduce the spin-glass phenomenology. In particular, dynamical rejuvenation effects 
in the absence of temperature chaos are observed. This rejuvenation is ascribed to the local 
dynamics, which is insensitive (due to the self similar nature of the potential) to the partic- 
ular valley that has been reached during aging at a higher temperature. Still, the separation 
of time scales/length scales with temperatures is much weaker that in experimental spin- 
glasses partly due to the rather modest time scales investigated in the present study. 
Correspondingly, abrupt rejuvenation as the temperature is decreased and strict memory 
when the temperature is cycled cannot be achieved. Rejuvenation and memory are present 
in embryo. 
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In a word, the Sinai model is a smeared out version of the multi-level trap model studied 
1^, 21 , where the sequence of critical temperatures is replaced by a gradual freezing in the 



details of the fractal landscape. One could have studied the effect of temperature shifts in a 
one dimensional potential with logarithmic correlations (rather than linear, as in the Sinai 
case). This model was studied in details in and was shown to exhibit a true transition 
temperature. We expect the resulting temperature shift phenomenology to be very close to 
that of the Random Energy Model, explored in ||2^. An extension of this model, that has 
an infinite sequence of phase transitions, will be presented in the near future |57| . 
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Appendix A: Calculation of the Entropy Fluctuations in the deepest well. 

In this appendix we compute the average entropy fluctua tions of the Sinai model in the 
large L limit. As we have already said in section IIIBl these entropy fluctuations are 
constant in this limit. For this computation we only need to consider the contribution to 
the partition function arising from the deepest minima. In a recent paper 1 39 Monthus and 
Le Doussal compute the probability distribution of the partition function of the deepest well 
considering that the origin in energies is at the bottom of the potential. They compute P{Z) 
obtaining that the Laplace transform is 

' 1 2 

/2 (a y/t) 13a 

where P = 1/T and a^ is the variance of the random force. This result implies that we can 
work in terms of a dimcnsionless partition function z = Z/a'^ to obtain, 

1 

dze-'M^) = — -=- . (A-2) 

Since the energy of the absolute minimum is set to zero, the entropy fluctuations can be 
directly computed as, 

^{jy - (W)^) = In' Z - hTz" . (A-3) 

Now in terms of the dimcnsionless partition function we have that In Z = In z + In , thus 
when computing fluctuations of averages over z terms which depend on a (and thus on T) 
cancel to yield, 

'S(Tf - [SiT^^ z-\^ ^ j dzp{z) {\nzf - ij dzp{z)hiz\ , (A-4) 

which is a constant as we expected. This constant can be evaluated by computing the aver- 
ages of In z and In^ z. The starting point of the calculation is Derrida's integral representation 
of the logarithm. 



In. 



1 

vZ^Yiiaj - (e-* - e"*^) , (A-5) 

When we perform the average of e~*^ we obtain the Laplace Transform given in ( |A-2| ) thus 

1 / 1 



\nz^ lim / dt-\ ^rf^ - e-* . (A-6) 
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This integral can be split into two parts, the exponential integral 
°° 1 

dt -e * = -7 - Ingr + gr + -|- + ... lim gr ^ , (A-7) 
i 2 

where 7 is the Euler constant, and the part with the Bessel Function, 



t /2 

it) J^, V dt Io(t) J loiV^) 



= - Ing + 21n2 - 27 + {^{2) - l)f + ... (A-8) 



where we have used the Wronskian property of the Bessel functions 

-^ = Kit)Ut)-K^it)m . (A-9) 
Adding this two contributions we obtain that the average entropy reads. 



S'(T) = lnZ = -21na-21n2 + 7 . (A-10) 
To compute the average In^ z we have to evaluate the following integral. 



\n^z= lim / dt dt' — e-*-*'--^-^ '— ., 1 . (A-11) 

The contribution arisi ng fro m th e first three terms in the integral can be easily evaluated 
from expressions (A-6) and (|A-7|) to obtain: 



dt I dt' ^ — = J" + "fin qp + hi q hip + 0{q,p) 



00 ^00 g— t ^ 



''I ''TV1!W) ^ (^ + l-'^)(27-21n2 + lnp) 

r°° r°° e~*' 1 

J dtj dt' — -^^^ = (7 + lnp)(27-21n2 + lng) . (A-12) 
The last term is somewhat more tricky and can be evahiated as follows, 



dtj dt'-^^2 rdt'-^^+2 rdt''-,^^ .(A-13) 

In the limit p — > the integral can be Taylor expanded yielding, 

dt'-^^^Tdt'-^+Oip) (A-14) 

where the t erms of order p do not contribute when we set p to 0. Hence we can evaluate 
expression ( A-13 ) at p = to obtain, 

/ dt i^^^ = -ln2g-2 ln21ng + 27 lng + 2C C = 0.2415... (A-15) 

Adding all the contributions and taking the double limit g — > and p — s- we obtain the 
following result for In^ z, 



ln2z = 4C-37 + 47 ln2 , (A-16) 
which finally yields the entropy fiuctuation, 

K^S^-S'^=A (C- (7- In 2)^) = 0.912784.... (A-17) 
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Appendix B: High temperature expansion of the distance shift 



We want to compute the high temperature behavior of the average distance shift between 
systems at two different temperatures di^2 a-s defined in eq.(^, 



((^>i - (^)2)2 (B-1) 
fn the high temperature hmit /? — > a simple Taylor expansion yields, 

(.) = ^ (e - (E -^^^ - E E t) ) (^-2) 

where Vi is th e ran dom potential at position Xi = i/L. 
Therefore for ( B-1 ) we obtain, 



{.{x)i-{x)2? = A/3';^(^E^^^^"E^*E|) 



, (£^-l)(L^ + 10L + ll)) 
120 1? ' ^^"^^ 



where in the last equality we have used the relation Vi Vj — cr^ min{i,j} which holds for 
Gaussian distributed forces with zero mean. For large L we recover expression (pf). 



i{x),-{x)2r = —iAprL (B-4) 
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